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Abstract. Following Nye and Berry's analogy with crystal dislocations, an approach 
to the Burgers vector of a wave dislocation (phase singularity, optical vortex) is 
proposed. It is defined to be a regularized phase gradient evaluated at the phase 
singularity, and is computed explicitly. The screw component of this vector is naturally 
related to the helicoidal twisting of wavefronts along a vortex line, and is related to the 
helicity of the phase gradient. The edge component is related to the nearby current 
flow (defined by the phase gradient) perpendicular to the vortex, and the distribution 
of this component is found numerically for random two-dimensional monochromatic 
waves. 

PACS numbers: 03.65.Vf,42.25.Hz,61.72.-y,47.32.~y 
1. Introduction 

In their study of phase singularities as a general phenomenon in wave interference and 
diffraction, Nye and Berry called them wave dislocations [1]. Other names used for these 
objects include optical vortices (as they are circulations of the flow associated with the 
waves), and nodal points in two dimensions, and nodal lines in three, as they are also 
the zero loci of the wave intensity. 

The term 'dislocation' is more established for defects in crystals, that is, 
imperfections in a regular periodic lattice geometry, which also occurs at points in two 
dimensions and on lines in three, and were originally proposed as a general feature of 
real crystal lattices in the 1950s [21 [3] 0. Since then, the study of geometric topological 
defects has played a key role in condensed matter physics and beyond (e.g. O [6]). 

These defects are topological because their presence signifies a breakdown in a 
key geometric descriptor of a spatially extended field (the order parameter), and a net 
change in that parameter in the vicinity of the defect, for instance on taking a circuit 
around it. For phase singularities, this is the strength: the phase change, quantized in 
units of 2tt, around the node. For nodal points in two dimensions, this gives a signed 
integer (topological charge), whereas for lines in three dimensions, the sign endows the 
line with an orientation (topological current): a natural direction for the tangent at 
each point along the nodal line. 

I It is interesting to note that around the same period, special examples of optical vortices were first 
noticed and studied [Ij. 
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Figure 1. Wave and crystal dislocations, (a) edge and (b) screw dislocations in 
crystal lattices, with arrows showing the Burgers vector directions, (c) and (d) show 
wavefronts (zero contour surfaces of real and imaginary parts) of the edge (c) and 
screw (d) dislocated waves described by Eq. , with propagation directions k given 
by the arrows. 



A dislocation in a crystal lattice, on the other hand, are defined by its Burgers 
vector |2j: the vector by which a circuit around the dislocation would fail to close in a 
corresponding circuit in a perfect lattice. It is thus a vector with integer components, the 
two most distinctive being an edge dislocation, whose Burgers vector is perpendicular to 
the direction of the dislocation line (all dislocations in 2-dimensional lattices are of this 
type), and a screw dislocation, whose Burgers vector is parallel to the line's direction. 
Examples of these two kinds are shown in Figure [D^a) and (b). 

Crystal and wave dislocations are described by different mathematical objects since 
the fields of which they are defects have different physical nature - a wave field, denoted 
here by ^Z', is a smooth, complex function of position, and a crystal is a discrete lattice 
generated by specific lattice vectors. Thus, does it make sense to define a Burgers vector 
for a wave dislocation? 

Nye and Berry's choice of nomenclature was motivated by specific, simple examples 
where dislocated waves have similarities to their crystal counterparts. Specifically, they 
considered a time-independent wave function of the form 



with a node where x = y = (giving a point at the origin in two dimensions, and a line 
along the 2;-axis in three). This is modulated by a plane wave with wavevector k. With 
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k in the x, y-p\a.ne (say the y-direction), the pattern of constant phase hnes resembles a 
crystal edge dislocation. In three dimensions, if k is parallel to z, the surfaces of constant 
phase are helicoids and resemble a screw. Examples of these surfaces of constant phase 
in three dimensions are shown in Figured] (c,d). Helical screw dislocations occur on 
the axis of optical beams carrying orbital angular momentum, such as Laguerre-Gauss 
modes [3 [8]. Mixed edge-screw dislocations occur as intermediates. 

However, general wavefields, typified by superpositions of plane waves propagating 
in different directions, do not in general have an overall carrier wavevector, although 
they do have a complicated configuration of nodal points or lines. Nevertheless, at 
points where the phase is not singular, it is natural to define the local wavevector as 
the local phase gradient [9], [TOl [11]. In quantum waves, this definition is related to the 
concept of weak measurements (i.e. local rather than global averages) [12]. Here, I will 
consider a general definition for the Burgers vector of a wave dislocation, based on a 
regularized phase gradient at the phase singularity. 

This work extends my previous approach [13] to the local phase structure of a 
wave dislocation. There, I considered several different measures for the rate of twist 
of the phase helicoids around a general dislocation line. This twist geometry is quite 
complicated, since the intensity contours in the plane transverse to a dislocation line 
are ellipses, with associated nonuniformity in phase [T^ [T5| [16] . and generically this 
anisotropy ellipse rotates - 'twirls' - along the singular line, resulting in the different 
phase helicoids twisting at different rates. This was used to define an average rate 
of rotation, the 'ellipse-defined twist', of the phase core along a singularity line ( [T3] 



where •' = t ■ V», and t is the tangent vector of the dislocation line (whose direction 
is provided by the topological orientation), parallel to ImVip* x V'?/' [I8j. As a measure 
of screwness, this quantity has certain attractive properties, although its derivation was 
somewhat oblique: it is a pure derivative of a certain angle, implying that it integrates 
to an integer times 2tt around a vortex loop (the screw number of the dislocation), 
and, in random waves, its fluctuations were smaller than the other measures considered. 
I will show here that Tweii arises naturally as the tangential screw component of the 
wave dislocation's Burgers vector, and find the corresponding transverse part. Other 
measures considered in Ref. [13], such as the azimuth-averaged twist Tw^ (the average 
rate of rotation of all phase surfaces along the singularity), and the screwness a (the rate 
of change of phase for fixed azimuth), appear as the tangential components of alternative 
definitions of the Burgers vector. 

Nye proposed [17] a similar definition of the Burgers vector for dislocations in time- 
dependent waves (based on a generalization of tpk), although it is difficult to generalize 
his definition to time-independent waves. In particular, his definition is based around a 
specific coordinate system attached to each point on the singularity, based on the tangent 
to the singularity and its perpendicular direction of motion. Nye's Burgers vector is 



Eq. (2.8)) 



TWell 
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then the hmit of the phase gradient on approaching the singularity from these specific 
orthogonal directions; when the dislocation is not moving, these directions cannot be 
defined. The Burgers vector defined in this paper is based on a local average of the phase 
gradient at points in a neighbourhood perpendicular to the dislocation line; the sense 
of circulation around the phase vortex cancels out in this average, leaving the residual 
phase gradient which defines the Burgers vector. The waves considered here will be 
time independent and monochromatic (although these assumptions are not completely 
necessary) . 

The structure of this paper is as follows. In the following section, basic notions 
associated with phase gradient and phase singularities will be described, along with 
the quantum-style Dirac notation required to write the Burgers vector. The derivation 
and basic description of the Burgers vector is given in Section [HI culminating in its 
definition in Eq. (13.31) . The details of the calculation here is given in the Appendix. 
Readers interested only in the consequences of the definition may omit these sections. 
The direction of the Burgers vector in the direction of the dislocation tangent - the screw 
component - determines the helicoidal twisting of the phase surfaces along the line, and 
this is related to the notion of helicity in fluid flows in Section HI The remainder, a 
vector perpendicular to the singularity line, defines the edge component, is discussed in 
Section [51 In particular, the statistical distribution of the edge component in isotropic 
2-dimensional monochromatic waves are explored numerically. 



2. Local wave dislocation geometry 

The following discussion applies to a general complex scalar wavefunction 

?/> = pexp(ix), (2.1) 

with intensity p^, and phase x- Explicit dependence on position (either in two or three 
dimensions) will usually be omitted. Wave dislocations are loci where ip = 0; on these 
points (2D) or lines (3D), the intensity vanishes quadratically, and x is undefined. As 
vortices, the dislocations are also circulations of the current associated with the wave, 

j = lmij*Vip. (2.2) 

The wave Burgers vector, defined in the next section, will be expressed in terms of 
various derivatives of the field close to the dislocation, which are described in this 
section. Similar discussions of this geometry may also be found in Refs. [H [TSj [TBI [19] . 

The phase is a well-defined function at typical points off the nodes; its gradient 
may be written 

Vx = ImVlog-?/' = Im-^^ = (2.3) 

Vx is the local wave momentum at the point evaluated, i.e. the local wavevector p| [TU]. 
For waves described by a vector field (such as the electric vector for polarized light), the 
natural definition for local propagation direction [20] can be written as a local wavevector 
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in an analogous way; these quantities for non-singular points are explored in more detail 
by Berry |Tl]. When ip is zero, however, Vx is singular, since if) and p are zero - this 
definition for a local wavevector fails. 

As described in the introduction, the topological strength s of a phase singularity 
is determined by a line integral of phase on a closed circuit C around the singularity. 
Using Stokes' theorem, 

s = ^ i dr ■ Vx = / dV ■ (V X Vx), (2.4) 
2vr Jc 2,7!: Jg 

where S is the surface or area enclosed by C, and s is the integer representing the net 

number of circuits of phase executed with respect to C; it will be assumed throughout 

this paper that the singularity has unit strength, i.e. s = ±1. Since Vx is a well defined 

gradient off the singularity, its curl vanishes there, and the only contribution to the 

integral in fl2.4p can be a ^-function at the node, 

Vx Vx = vr52(^)(Vx j). (2.5) 

(This may be more rigorously justified by regularizing Vx at the singularity.) Equation 
(12. 5p shows that the wave's vorticity - the curl of the fiow velocity Vx ^ is concentrated 
on the vortex lines with weighting 

i7 = X j = yV^* X V^. (2.6) 

O alone is often referred to as the vorticity, although as the curl of Vx (rather than j) 
it is only nonzero on the nodal lines. In three dimensions, 17 points in the vortex line 
tangent direction t in the direction of topological current; in two (denoted by x,y), its 
scalar sign (17 dotted with the unit 2;- vector) defines the topological charge sign. 

It is convenient to choose local cartesian coordinates defined by this direction. With 
the origin on the singular line, with tangent in the z-direction (parallel or antiparallel 
to the topological current), the 2-dimensional vorticity may be written 

17 ■ i = ^{r.^y - ry^.) = im ^vr ■ ^3 ■ v^. (2.7) 

Acting between the gradient vectors in the final equality is the Pauli matrix cr^ = 

-i 

1 

that the labels are permuted from their usual place, as they are here applied to vectors 
in a cartesian basis, rather than a polar basis), for j = 0, 1, 2, 3, give 

Sj = V^* ■ aj ■ Vip. (2.8) 

These Stokes-like parameters describe the local anisotropy of the gradient vector Vijj 
[T3t [T5l \T6\ [19] . Thus Sq = jVV'p, and S^, = 217 ■ i, whose sign gives the singularity 
strength. The orientation angle of the anisotropy ellipse defined by Vtp is given by 
|arg(S'i + 15*2). 

The anisotropy ellipse describes the elliptical contours of low intensity, and related 
squeezing of phase contours around the vortex; since the fiow of current j = p^Vx is a 
perfect circulation, these two effects cancel [TH [H]. The major axis of the Vip ellipse 
corresponds to the minor axis of the local intensity ellipse, and vice versa. 



. Similar inner products with respect to the other Pauli matrices (noting 
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Vy = — Im I — I -|- ini 



0{R), (3.1] 



3. The Burgers vector as a regularized phase gradient 

In this section, the Burgers vector will be defined following the definitions and notation 
of the previous section. With coordinates chosen so the dislocation is at the origin 
(with direction in the z-direction for three-dimensional waves, and the z = plane is 
considered), the phase gradient Vx is singular as the cylindrical radius R = a/x^T^ 
approaches 0. Around the dislocation ip = 0, at which the field gradient is V^ipo, the 
phase gradient Vx = Im.'Vtp/tp may be expanded in powers of R, 

' ^ u Vipo u ^ u ' 
u ■ VipQ 2~(n- V^o)^ 
where u denotes the unit radius vector u = {x,y,0)/R = (cos 0, sin0, 0) and ^ 
is the hessian matrix of second derivatives, "^ij = didjip, i,j = x,y,z. The first 
term of (13 .ip is proportional to R~^, accounting for the singularity. The numerator 
of this term proportional to the direction field for the current near the singularity, 
j = Rim {{Vip ■ u)'Vip*) + 0{R^), known to be a perfect circulation [18], [H]. 

The second term is independent of i?, although dependent on unit direction u. Its 
two components (within [ ] brackets) arise from the next-to-leading order terms in the 
expansion of '^ip/ip, and depends linearly ion the second derivatives ^ . Only the first 
of these has a component in the 2;-direction, as the other is proportional to VV^o, which 
is by definition orthogonal to z. The expression in [ ] accounts for Vx at the singularity 
(up to terms of order i?), aside from the singular circulation in R^^. 

Evaluating this -R-independent term for the dislocated waves (11. ip indeed gives 
the carrier wavevector fc; this in fact is also the case for the more general form 
Ip = {{X + iY) ■ r) exp(ifc • r), for X, Y any pair of linearly independent vectors in the 
X, y-plane. For these simple cases, the local wavevector is simply the non-circulatory part 
of Vx in the limit R 0, and is therefore a natural definition of a local wavevector. 
For more general waves, however, this is dependent on u, and so this limit must be 
approached democratically from all directions. 

The Burgers vector b for a wave dislocation in a typical field, therefore, is defined 
to be the integral of the phase gradient (13.11) over a small circular disk in the xy-p\ane 
centred on the singularity, in the limit of the disk radius d — > 0, 

6 = lim 4^ / dV Vx = 7^ / del) Vx- (3.2) 

The singular circulatory term cancels (since if) is averaged before R —>■ 0), so 6 is 
determined purely by the term in square brackets [ ] in (13. ip . This definition is more 
general, but similar, to Nye's definition of a Burgers vector in time-dependent waves 

These integrals, involving the ratios of linear and quadratic forms in u, are 
straightforward to integrate using complex contour integration; this is described 



in Appendix A The singular term in i? integrates to zero, constituting the 
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regularization. Explicitly, the integral of the terms in [ ] in fl3.ip gives 



Im 



0) 



\-2 



(3.3) 



where is the isotropic circular vector with the same sense s as the singularity (from 



Bs = X + siy. (3.4) 

as defined here may be multiplied by an arbitrary complex number (e.g. normalizing 
it) without affecting b. The two terms of b arise directly from integrating the two terms 
in [] in (D. 

Equation (13. 3p for the Burgers vector 6, as a regularized phase gradient at the phase 
singularity, is the main result of this paper. Geometrically, 6 is a measure of the local 
curvature of phase lines or surfaces local to a dislocation; precisely how b as defined is 
related to this curvature is not clear. A short calculation of the ^-component 6^ gives 
—Tweii, the negative of the twist defined in (II. 2p . The significance and meaning of this 
will be discussed in section IH The edge component perpendicular to the vortex, in the 
xy-plane, will be discussed in section [51 

The choice of regularizing by integrating over a disk centred on the singularity is 
natural, but not unique. Another obvious choice is to integrate over the anisotropy 
ellipse. Effectively, this may be done by replacing u with Re (e~'^^V'?/'o) and integrating 
with respect to 6. This gives a phase- weighted Burgers vector similar to b, but with e* 
in (13. 3p replaced by Vip^; the third component of this vector, along the direction of the 
singularity line, is minus the screwness a defined in Ref. [13| Eq. (2.4). Averaging with 
respect to the intensity anisotropy ellipse, whose axes are orthogonal to those of VV^o; 
is achieved by using the vector Re {e~^^^ {—ipoyX + ipoxV))- The result is again the same 
as (13. 3p . but with with e* replaced by —iplyX + 'ipQxV^ ^"^^ ^^e component along the 
singularity line is the negative of the azimuth-averaged twist Tw^ of Ref. [13] Eq. (2.3). 
Each of these alternative definitions give the desired vector k for singularities of the 
form (II. ip . although they are different for more general waves. It appears most natural 
to define the Burgers vector weighting all spatial directions around the singularity 
equally, as above in (13. 2p . (13. 3p : in particular, the alternative definitions do not have 
the attractive property of being related to helicity, as described for b in the following 
section. 



4. The screw component 

The screw nature of a wave dislocation line in three dimensions, as the magnitude of 
the Burgers vector in the direction of the singularity line t ■ b = —Tweii, is a measure 
of the rate of rotation of the phase structure along the dislocation line. The — sign can 
be justified from the well-known fact that for a screw dislocation of the form (II. ip . the 
helicoid is left-handed when 17 ■ > 0, and right-handed otherwise; Twew was defined 
in [13] to be positive for a right-handed screw. 
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This quantity, integrated over a simply-connected volume V whose surface is not 
pierced by any singularities, it is reminiscent of the helicity if of a fluid flow with velocity 
field V, 

H = / d^rv-V XV. (4.1) 
Jv 

Helicity is a topological invariant: provided the topology of the vortex lines does not 
change, H is constant. With v = Vx, the integrand is 27T6'^{tp)f2 ■ b, so the helicity of a 
scalar optical field is the integral of the screw component of the Burgers vector —Tweii 
along all of the closed loops contained in the volume V. 

Around each such loop, the phase pattern must rotate an integer number of times, 
which may be thought of as the topological linking number of the edges of any local 
phase ribbon centred on the singularity [131 122]. For a planar loop, this linking 
number is simply the integral of i7 ■ 6 around the loop, and if the loop is non-planar, it 
is the twist Tw for the Calugareanu- White- Fuller relation Lk = Tw + Wr (with linking 
Lk and writhe Wr) [221 [23] . The helicity H is also related to this theorem in more 
general fluid mechanics [23]. 

The significance of the previous discussion is that a closed dislocation loop with 
nonzero helicity - for which there is a net integer number of screw rotations of the phase 
around the loop - must be threaded by other dislocation lines [251 126] . Knotted (self- 
linked) dislocation lines have a helicity linking number given by their minimum crossing 
number . This result is a consequence of the continuity of the complex scalar field in 
three dimensions: the surfaces of constant phase must fill all space, only crossing each 
other on the singularity lines, and all phase surfaces must cross there. 

A simple example of a wave (solving the Helmholtz equation) with such a threaded 
loop is [251 [13] 

■^threaded = -R™ cxp (ims0) (A;(i?^ - a^) + 2i(l + m)z) exp{ikz), (4.2) 

in cylindrical coordinates {R,(f),z), and integer m > 0. This field has a straight 
dislocation line of strength m along the z-axis, with topological current parallel to 
sz, encircled by a nodal loop of radius a in the z = plane. This loop has both 
screw and edge character. The topology requires the integral of the screw component 
around the loop to be 27ism, and this was calculated in [13] to be —sm/a. The loop has 
edge nature since, lying in the xy-plane, it is perpendicular to the global propagation 
in z] its edge component is A; (1 — (1 + m)/2(l + m — ak^) in the 2;-direction. More 
complicated threaded and linked vortex lines have been seen to occur generically in 
computer experiments of random optical speckle fields [28], and knotted and linked, 
threaded vortices in holographically-controUed laser fields |29j . 

The general relation between vortex linking and knotting and helicity was described 
by Moffatt [30], and the discussion there for ^-concentrations of vorticity applies directly 
to the case of complex scalar wavefields. Helicity has previously been discussed along 
similar lines to here for wave dislocations and optical vortices [311 132] . 
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There are two simple cases of optical waves which are often described to contain 
'screw dislocations': paraxial beams (satisfying V^ipp = —2ikdzipp), such as Laguerre- 
Gauss modes |H]; such waves represent an approximation to paraxially small deviation 
from an exact plane wave travelling in the z-direction exp(iA;2;)^p. In the paraxial 
approximation, k is assumed to be infinitely larger than the largest transverse Fourier 
component in ijjp, so all the dislocation lines are approximately parallel to the 2;-direction. 
Nonparaxial plane waves may also have a constant ^-component exp{ikzz)ilj^fi^, and these 
are the nondiffracting beams whose intensity is unchanged on propagation in z, such as 
the Bessel beams [33]. It is easy to see in both these cases that the screw component 
of the wavevector, as defined in section [3], is the constant magnitude of the component 
of the 2;-component of the wavevector: k^ for nondiffracting beams, and for paraxial 
waves, approximately k. 



5. The edge component 

The edge component of the Burgers vector is a vector defined for all optical vortices 
(possibly zero for pure screw lines). For two-dimensional fields, such as the transverse 
plane of a (paraxial or nondiffracting) beam, this vector lies in the transverse plane, and 
apparently does not correspond to any previously described quantity associated with 
optical phase singularities. Unless stated otherwise, 'the Burgers vector' in this section 
refers to the edge component of the Burgers vector. 

The vector is the direction of the carrier wave local to the vortex line. For two 
simple examples with a global carrier wave exp(i%), namely the edge dislocated wave 
ipky of (11. ip . and the exact two-dimensional solution of the Helmholtz equation [1] 
^pair = {k^x"^ — a'^) + iy) exp{iky), are shown in Fig.O The latter case has a pair of nodes 
with s = ±1 at (±a, 0), annihilating when a —>■ 0. Analytically, the Burgers vector of 
these two dislocations is (0, A;(l — l/(H-2a/c)^)). Although in the limit a — > 0, this is finite. 



the strength s = at the annihilation event, and the assumptions in Appendix A do not 
hold. More general forms of annihilating dislocation points, as considered for instance in 
Ref. [Hlj, have Burgers vectors which are the sum of any overall wavevector and a term 
related to the vortex anisotropy ellipses. In three dimensions, the axisymmetric analogue 
of ^/'pair is an edge dislocation loop perpendicular to the propagation direction, and 
may be realised by setting m = in (14. 2p . for which the (purely edge) Burgers vector 
has length k{l - l/2/(-l + ak)^). 

In the two-dimensional edge plane perpendicular to the vortex line, there is a 
saddle point of the phase and flow of j (i.e. x finite, V x j < 0), the Burgers 
vector is perpendicular to the displacement vector between the node and the saddle. 
This phenomenon is straightforward to explain in the case of a dislocated field with a 
global wavevector, such as (II. ip : in the vicinity of the singularity {R small), the main 
contributions to the vector field Vx are the first two terms in (13. ip : the circulation term, 
proportional to R~^, and the Burgers vector term, constant when the wavevector is 
global. Very close to the vortex, the circulating term dominates Vx, but as R increases. 
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(a) (b) 




Figure 2. Edge dislocated waves in the plane, (a) -^fc^ from Eq. (|l.ip . and (b) V'pair- 
Colours denote phases labelled by hue, and the white curves show the streamlines of 
current and Vx, which tend to follow the global wavevector in the y-direction. The 
Burgers vectors, also in the y-direction, are the black arrows, represented four times 
their true length. Saddle points are plotted in grey. 




Figure 3. Four square wavelengths of an isotropic random wavefield, represented by 
the superposition of 100 plane waves with random directions, amplitudes and phases. 
The quantities represented are the same as in Fig. [2j 

the constant Burgers vector term has more effect. At the point in this neighbourhood 
when the two terms have equal magnitude and opposite direction, in direction s{—by, h^) 
with respect to b = (5^., hy) (assuming pure circulation), Vx = 0, i.e. at the saddle point. 
Although the Burgers vector term is not constant for more general fields, it is natural 
to assume that phase saddle points occurring near dislocations [3ll [35] are related to 
the direction of the Burgers vector in similar ways. 

A representation of a more complicated planar wavefield, which does not have a 
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Figure 4. Probability of the Burgers vector magnitude b for two-dimensional isotropic 
monochromatic fields, (a) Probability density P{b) represented by histogram (with 
horizontal axis units of k). (b) Cumulative probability C(&), with insets showing log- 
log plots of C{b) for small b (lower left) and 1 — C{b) for larger b (upper right). The 
linear fits correspond to 35. 9&^'^ for small b, and .0016&~'*'^^ for b > 0.8. b was calculated 
here from over 500 000 different dislocation points from large areas of several different 
realizations of superpositions of 200 random waves. 



globally defined wavevector, is shown in Fig. [3l The field depicted is a superposition 
of 100 plane waves with constant wavenumber k, propagating in uniformly random 
directions with complex circular gaussian random amplitudes. As is well known, there 
are many dislocations and saddle points (each with density k'^/Air [18j). Also plotted are 
the stream lines of current and Vx, and the (edge) Burgers vectors at the dislocations. 
Most of the dislocations in the figure have a nearby saddle point, situated to the left 
of the dislocation with respect to the Burgers vector, as suggested by the preceding 
argument. As can be seen on the upper right hand side of the figure, when there is a 
strong sense of direction to the current near isolated dislocations, the Burgers vector on 
the dislocation has a greater magnitude, and tends to be in this direction. When the 
nearby current does not have a strong consensus direction, as in the lower right half of 
the figure, the Burgers vectors are smaller. It is difficult to provide a more quantitative 
description of these phenomena in plane wave superpositions due to the form of b, which 
is very complicated even in the simple case of superpositions of three plane waves. 

A natural question to ask is what the statistical distribution of the Burgers vector b 
is for random waves of the type shown in Fig. [3l given the wealth of knowledge of other 
statistical features of wave dislocations (e.g. [HI |36l H]). In particular, the probability 
density function of the screw component of the Burgers vector for isotropic random 
three-dimensional waves was computed in [13] (Eq. (4.1)) to be a generalized Cauchy- 
type distribution. However, the complicated form of the edge part of b has so far proved 
intractable to statistical methods. The results of numerical computation of the length 
^ — \/^x + of the Burgers vector for isotropic 2-dimensional monochromatic waves (or 
random nondiffracting waves in three dimensions) is shown in Fig. H] as a histogram for 
the probability density function P(6) in (a), and the cumulative probability distribution 
C(6) = Jq P{b')db' in (b). From this numerical data, it is possible to estimate the root 
mean square Burgers vector length (62)i/2 as 0.22; clearly, the expected Burgers vector 
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length is somewhat smaller than k. Linear fits (shown in the insets to (b)) suggest that 
P{h) oc h for 6^0, with power-law decay asymptotically, P{h) ~ 0.0076~^'^ for large h. 
This suggests that there is a small but nonvanishing fraction of dislocations in this model 
for which the regularized local phase gradient b is larger than the global wavenumber 
k : numerically, this fraction is approximately 0.0016. This contrasts with the fraction of 
nonsingular points whose local phase gradient magnitude exceeds /c, which is 1/3 [10]. 

6. Discussion 

All generic nodes/phase singularities in wave fields are vortices of energy, regardless of 
their specific geometry or orientation with respect to propagation direction. However, 
other properties, such as the well-known helical wavefronts of optical screw dislocations, 
depend on the relationship of the zero with the global and local phase gradient, and 
this is summarised mathematically by the local phase gradient, here equated with the 
Burgers vector of the wave dislocation, in the spirit of Nye and Berry's approach. 

This local phase gradient, represented in equations (13.31) . depends both the first 
and second derivatives of the field at the node. As such, it is a second order property of 
phase singularities, whereas quantities such as the vorticity (12.71) and other Stokes-like 
anisotropy parameters (12. 8p only involve first derivatives. Nevertheless, the Burgers 
vector plays an important role: topologically, in that its screw component, determining 
the local helicoidal twisting of the wavefronts, is related to the helicity of the wavefield; 
and physically, in that the edge component is a measure of the strength of current flow 
beyond the immediate vortical circulation. 

This phase gradient is not the only quantity going beyond the immediate 
neighbourhood of optical vortices; for instance, the current streamlines themselves have 
interesting spiralling properties beyond first order expansions near vortices [371 138] . 
However, the spiralling defined in Ref. [37] depends on third as well as second derivatives 
of the field close to the vortex, and so is not immediately related to the quantities 
discussed here. 

Only basic properties of the Burgers vector have been discussed here, and further 
examination might yield further insights into optical singularity physics and the 
relationship between local and global energy flow near vortices. For instance, can the 
edge component of the Burgers vector twist along wave dislocation lines? If so, does 
this define a higher-order helicity for closed loops? 

The Burgers vector here defined might be related to the various phase gradient 
methods (e.g. [39l HQ]) proposed for the propagation dynamics of vortex points in 
paraxially propagating fields, which account for the transverse path of vortex points 
(effectively the tangent direction in three dimensions) by constructing a local phase 
gradient. Such a direction is either perpendicular ('glide') or parallel ('climb') to the 
direction of the Burgers vector (Nye and Berry [Ij studied such phenomena in time- 
dependent fields), and such a comparison with the crystallographic case may provide 
new understanding of the propagation of optical singularities. 
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Appendix A. Contour integration of the regularized phase gradient 

Integrals over < < 27r of functions such as 

= — Ah ■ (A-i) 

a cos + sm 
. ,s ccos0 + rfsin0 
a cos + sm 

, , , . c cos^ (p + d sin^ + 2e cos sin / . \ 

^ = 7 , , , . ,X2 ) A-3 

(acos0 + osm 0)^ 

with a,b ^ 0, may be evaluated using complex contour integration methods. 

In the following, g{(j)) will be integrated as an example. In doing so, it is useful to 
introduce the following function of complex variable Z, 

^ c{Z + 1/Z) - idjZ - 1/Z) ^ Z\c-id) + {c + id) 
^ ' ~ a{Z + l/Z)-ih{Z -l/Z) Z^{a-ib) + {a + iby ^ 

since G(e"'^) = 5'(0). Around the unit circle C in the complex Z-plane, dO = —idZ/Z, 
so 

"27r r r]7 

d4>g{4>) = -id> —G{Z). (A.5) 



'0 Jc Z 

G{Z) is a meromorphic function of Z, with simple poles at Z± = 
{a + ib) /(a — ib). These poles are inside the integration contour if the sign 

S = sign (|a + i6p — \a — i6p) = signlma*6 (A. 6) 

is positive; when g{(f)) corresponds to the -R-independent term in (13.21) . S = 
signlm'0*'0y = s, the singularity strength. In this paper, it is always assumed that 
this is nonzero; otherwise there are poles on the integration path in the complex plane, 
and higher-order regularization of the phase singularity is necessary. 

If S" = —1, then neither pole of G{Z) is enclosed by the contour, and the only 
contribution to the integral (lA.SP is the residue of the simple pole at the origin, giving 

"^'^ r-\-id 

d0 g{(P) = -i X 27riG'(0) = 2%—— (for S < 0). (A.7) 







a + ib 

When S = +1, the conjugate counterpart G{Z*) to G{Z) can be defined. 



^,^,,_ Z*^c + id) + {c-id) 

^ = Z*\a + ib) + ia-iby ^^-^^ 
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with g{(f)) = G{e~^'^), and on the unit circle in the Z'*-plane, d(f) = idZ*/Z*, and the 
integration contour for the integral corresponding to flA.5l) is clockwise. 

The simple poles of G{Z*), considered as a meromorphic function of Z*, are outside 
the unit circle, and do not contribute to the complex contour integral, which depends 
on the residue at the pole at Z* = in the usual way. Thus, for any g{(p), 

/'27r p _ j 9^ 

/ d0^(0) = 27r — , ^ = signlma*6. (A.9) 

Jo « - 

By similar arguments, the meromorphic functions F{Z), F{Z*) corresponding to 
have poles in the same places; however, the integrand has no pole at the origin, 
and so 

d0/(0) = O. (A.IO) 



The integration of h{(f)) follows a similar argument, and the corresponding 
meromorphic functions H{Z), H{Z*) with = H{e^'^) = HieT^'^) now have double 
poles at the same places as the corresponding functions H{Z), H{Z*). The arguments 
to integrate these are identical to those above, giving 

"27^ p _ ^ _ 2i S'p 

d<j)h{<f)) = 27r- — — , 5 = signlma*6. (A.ll) 

(a — iDo)^ 
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